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Abstract 

We construct new solutions of the vacuum Einstein field equations in four dimen- 
sions via a solution generating method utilizing the SL(2,i?) symmetry of the reduced 
Lagrangian. We apply the method to an accelerating version of the Zipoy-Voorhees 
solution and generate new solutions which we interpret to be the accelerating versions 
of the Zipoy-Voorhees generalisation of the Taub-NUT solution (with Lorentzian sig- 
nature) and the Zipoy-Voorhees generalisation of the Eguchi-Hanson solitons (with 
Euclidean signature). As an intermediary in the solution-generating process we obtain 
charged versions of the accelerated Zipoy-Voorhees-like families of solutions. Finally 
we present the accelerating version of the Taub-NUT solution and discuss its properties. 
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1 Introduction 



Ever since the formulation of General Relativity, exact solutions have played an integral part 
in our understanding of the nature of spacetime. For example, much of our understanding 
of black hole thermodynamics and inflation were possible only with the discovery of the 
Kerr-Newman and FRW solutions respectively. Given the importance of such exact solu- 
tions, there is a corresponding impetus to derive new solutions which upon analysis would 
yield further insight into our universe. Since Einstein's equations in their unadulterated 
form consist of a series of coupled non-linear differential equations, obtaining solutions by 
hand is intractable unless some kind of simplifying symmetry is imposed in the ansatz. This 
motivated the development of many ingenious and powerful strategies utilizing such simplifi- 
cations to derive solutions to Einstein's equations. Of particular importance for the present 
work is a special type of simplifying ansatz, the static axisymmetric Weyl-Papapetrou metric, 
which was first proposed by Weyl in pi 

dsl = -e-'^dt' + [e2^(c/p2 + dz^) + V] • (1) 

The metric is specified by the values of two functions ip and p, which are functions of the 
canonical Weyl variables p and z. Starting with a static vacuum axisymmetric solution of 
the Einstein field equations, consider its dimensional reduction along the timelike direction 
down to three dimensions. The reduced Lagrangian can be written as: 

£3 = e/? - ie(9^)^ (2) 

where we denote e = ^/g and R is computed with the three-dimensional metric: 

ds^ = e^^'idp' + dz^)+ p'd^^. (3) 

The equation of motion for ip is then readily seen to be A^/^ = 0, where A is the Laplacian 
constructed using the three-dimensional metric (jH)). Now, the key observation is that Aip = 
e~^^A^=oV') where A^=o is the Laplacian computed for a flat three-dimensional Euclidean 
metric, which corresponds to setting /i = in Q. Therefore any solution ipip^ ^) of Laplace's 
equation in the fiat three dimensional space is automatically a valid solution of Laplace's 
equation in the curved background Once we know ip, the remaining function iJ,{p,z) is 
found by performing a simple line-integral using the relations: 

d.p = ^-dp^pd^^P, = ^ [{dp^f - [d^^f] . (4) 

The Einstein's equations in 4-dimensions for a static axisymmetric background are now es- 
sentially reduced to finding a solution of Laplace's equation on fiat space. Due to the hnearity 
of the Laplace equation for tp, construction of multi-black hole versions is easily carried out. 
The Weyl formalism has been recently extended to higher dimensions by Emparan and Re- 
all j2] and the same line of thought can be used for the corresponding higher dimensional 
axisymmetric metrics. 

Given the simplifications introduced by the above axisymmetric ansatz, one now has two 
choices. One may either try to solve the differential equations directly, or, in more general 
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cases, try to further exploit the hidden symmetries of the dimensionally reduced Lagrangians 
and generate solutions using pre-existing solutions as seeds. 

Consider for instance the following 'scaling' symmetry of the field equations: given a 
vacuum static solution described by the pair of functions {ip,f^) then it is easily seen from 
dH that the pair {jip, will describe new vacuum static axisjTumetric solution of the 
field equations, where 7 is any real parameter. As an example, its application converts the 
Schwarzschild solution into the Zipoy-Voorhees solution This simple observation will be 
liberally utilised in the present work. 

There also exist transformations similar to the Ehlers- Harrison transformation for the 
Ernst formalism which map static vacuum solutions into stationary Einstein-Maxwell 

solutions. For the Weyl-Papapetrou ansatz, it has been long known that a transformation al- 
ready exists that brings a static, axisymmetric vacuum solution to a non-trivial class of static 
solutions in Einstein-Maxwell theory [7 . In particular, the Schwarzschild solution can be 
transformed into the Reissner- Nordstrom solution. In this paper, we will demonstrate a sim- 
pler alternative derivation of this transformation using a SL{2, R) symmetry of the reduced 
Lagrangian in three dimensions.^ However, unlike previous applications of this transforma- 
tion, we show that combining this transformation with the above scaling symmetry, we are 
able to generate new solutions. Our seed metric will be an accelerating version of the Zipoy- 
Voorhees solution Pj. In particular, we obtain new vacuum stationary axisymmetric metrics 
which we interpret as describing the accelerating Zipoy-Vorhees-like family of Taub-NUT 
solutions (with Lorentzian signature) and Eguchi-Hanson instantons (with Euclidean signa- 
ture). Much like the original Zipoy-Voorhees solution [3^, such metrics are parameterized by 
a real number a. For a = 1 we recover the accelerating Taub-NUT/Eguchi-Hanson solitons 
and, for higher positive integer values of a, they can be interpreted as the 'superposition' of 
accelerating a NUT-charged objects/sohtons. 

The structure of this paper is as follows. In section |21 we describe our solution generating 
technique, which maps a static axisymmetric solution in vacuum to a new stationary vacuum 
solution of Einstein's gravity in four dimensions. In section |21 we apply the transformation 
technique on accelerating solutions, namely on the C-metric and on the accelerating Zipoy- 
Voorhees metric and we consider more closely the properties of the generated accelerating 
Taub-NUT solution in section |3] Finally, in section El we summarise the main results of 
this paper and discuss potential avenues for further research. In Appendix A we cast the 
accelerating Taub-NUT metric in Weyl-Papapetrou form. 



2 Weyl's charging method: the SL{2^ R) approach 

We start with Einstein-Maxwell theory in four dimensions described by the Lagrangian: 

A = ei?-^eF2), (5) 

where R is the Ricci scalar computed with the 4-dimensional metric, F(2) = (iv4(i) is the 
electromagnetic field strength which only has an electric component = xdt ^ind we 

similar idea has been considered in however the details of the charging transformation differ here. 
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denote e = \/—g- Consider the dimensional reduction of the four- dimensional Lagrangian 
(0) to three dimensions on a timelike coordinate using the static Kaluza-Klein ansatz: 

dsl = -e-^de + e'^dsl (6) 

The reduced Lagrangian in three dimensions is then^ 

C, = eR-\e{d(t>f + \ee^{dxf. (7) 
Notice now that if one defines the matrix: 

M = 



y ± -± , y'^ ± 

— e 2 + 



then the three-dimensional Lagrangian can be cast into the following form: 

£3 = ei? + eTr [dM~^dM\ , (9) 

The reduced Lagrangian is then manifestly invariant under general SL[2, R) transformations 
if one considers the following transformation laws for the three-dimensional fields: 

g^.u ^ g^.u, M ^ Q^'MQ, n=(''^], ad-bc = l. (10) 



d 

Starting now with a static axisymmetric vacuum solution described by the metric: 

ds^ = -e-'^df + e'^dsl, (11) 

then performing the dimensional reduction on the timelike direction down to three dimen- 
sions and applying a general SL{2, R) transformation parameterized as above, we obtain 
a static axisymmetric electrically charged solution of Einstein-Maxwell field equations, de- 
scribed by the fields: 

ds"^ = -e-'^df + e^dsl, = xdt, 

where, in terms of the parameters appearing in Q, the new constants 6 and C can be expressed 
as 6 = /a? and C = l/(2ac). Note that in the limit in which Q = I2, i.e. c — * and a = 1, 
we have 5 — > simultaneously with C —>■ ±00 such that the product C^6 —>■ 1/4 remains 
constant. 

As an example of this charging technique, let us generate the Reissner- Nordstrom solution 
starting from the Schwarzschild metric: 



rfs^ = - ( 1 - dt' + + r\de' + sin^ edif^). (13) 



Here e ~ vM ^^'^ ^ computed with the 3-dimensional metric c?s§ 
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The final solution can be written in the form: 

{{l-5)r + 25mf AC^5r{r-2m) "^"^ + ^CH ^"^^ ^^'^ ^"^^ 

2C{{l + 6)r-26m) 

^(^^ - {l-5)r + 25m ^^^^ 

For generic values 5 7^ 1 we can easily perform a redefinition of the radial coordinate, together 
with an appropriate constant scaling of the timelike coordinate and cast the solution into 
the usual Reissner- Nordstrom form. The electric charge is Q = m/C, while the mass of the 
solution is M = {1 + 6)Q / {2\/6) . Note that 5 = 1 is a special case as it leads to the Bertotti- 
Robinson metric ^U] and it therefore describes the extremally charged Reissner-Nordstrom 
solution for which M = Q. 

At this point let us note that this method is not restricted only to metrics with axisym- 
metric symmetry; it can be extended to any general static vacuum solution of Einstein's field 
equations. Also, one can easily consider a similar method to generate magnetically charged 
static solutions out of axial vacuum metrics. 

Once we obtained a static electrically charged solution in Einstein-Maxwell theory, the 
next step in our solution-generating technique will be to perform a dualisation of the electro- 
magnetic potential and find the corresponding magnetically charged solutions. In our case 
it turns out that it is easier to compute the dual electromagnetic potential in the reduced 
three-dimensional theory. In this case we start with the three-dimensional Lagrangian ((7j) 
and dualise the scalar field to obtain a magnetic 2-form field strength F(2). Following the 
usual duahsation procedure, we add a term dx A F(2) to the action and solve the equations 
of motion for the scalar field x- Replacing the result in the action we finally express the 
Lagrangian in terms of the dual field as:^ 

£3 = ei?- ^6(90)2 -iee-^F(|), (15) 
where the components of the two-form field strength are computed using the formula: 

Fa/3 = ee'^ean^d^x. (16) 

where eQ,/3^ is the Levi-Civita symbol. After lifting the solution back to four- dimensions we 
obtain a magnetically charged static solution of the Einstein-Maxwell field equations. 

We are now ready to perform the last step in our solution-generating method, namely to 
map the magnetic solution to a vacuum axisymmetric stationary solution of Einstein's field 
equations in four dimensions. This actually involves two steps: we first map the magnetic 
solution of the Einstein-Maxwell theory to a solution of the Einstein-Maxwell-Dilaton (EMD) 
theory with a specific value of the dilaton coupling, namely the one corresponding to the 
Kaluza-Klein theory, i.e. a = —^/S. To do this we shall employ the general results derived 
in ^T] (see also [12j for a geometrical derivation of the respective mapping). Starting with 
a magnetostatic solution: 

dsl = -e-^df + e't' [e^^{dp^ + dz^) + p^d^"^] , = A^d^, (17) 



^Note that we perform the duahsation using a three dimensional Euclidean metric. 



5 



the corresponding solution of the EMD system is: 



dsl = -e-Ut^ + ei (e^^) ^ (dp^ + dz") + p^d^^ , 



(18) 



This is none other than the dimensional reduction of a vacuum five-dimensional metric using 
the ansatz: 



In our case it turns out that the 5-dimensional metric is simply the trivial product of a 
4-dimensional Euclidean metric with a time direction. Since the 5-dimensional metric solves 
the vacuum Einstein equations it is manifest that the 4-dimensional Euclidean metric will 
be Ricci fiat, i.e. it solves the vacuum Einstein equations in four dimensions. Therefore, our 
final result is expressed in the form: 



We note that even if the charging method does not require any other symmetry beyond 
the static condition, this second step in our solution-generating technique can be applied 
only to stationary axisymmetric metrics that can be cast into the Weyl-Papapetrou form. 
However, this is not really a very stringent constraint as most of the physically interesting 
solutions can be cast in the Weyl-Papapetrou form. To understand the effects of this last 
step in our solution-generating method one could take for instance the magnetically charged 
four- dimensional Reissner- Nordstrom and map it to an Euclidean vacuum metric as in (|2Up. 
However, given the presence of the square roots appearing in the factors ez it is easily 
seen that we obtain an axisymmetric NUT-charged solution with naked singularities, whose 
physical interpretation is obscure at this point. On the other hand, since we restrict ourselves 
to axisymmetric metrics that can be cast in Weyl form, it turns out that before we apply 
the charging procedure one can use the scaling symmetry discussed in introduction to scale 
the dilaton in the initial seed in such a way to cancel the awkward effect of the square-roots 
in the final expression of the metric. 

3 Accelerating Zipoy-Voorhees-like families of solutions 

In four dimensions, a particularly interesting class of solutions that generalise the Schwarzschild 
black-hole is the so-called C-metric. The static part of this metric was found by Levi-Civita 
almost one century ago (see [H]), however, its physical interpretation was clarified only af- 
ter Kinnersley and Walker's work decades later By performing appropriate coordinate 
definitions, they found that this metric describes a pair of causally disconnected black holes 
uniformly accelerating in opposite directions. The cause of the acceleration is understood 
in terms of nodal/conical singularities along the axis that connects the two black holes and 




(19) 



dsl = e"^ {dz + A^/2d^f + e 2 (e^^) ^ {dp^ + dz^) + p^d^"^ . 



(20) 
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these singularities are interpreted as strings/struts pulling or pushing the black holes apart. 
A more general class of electrovacuum spacetimes that includes and considerably generalises 
the C-metric was found by Plebahski and Demiahski jll4j (see also in 15] the general Type D 
family of solutions in EM theory). Recent analyses of this class of solutions have been per- 
formed in El , where a new exact solution describing a pair of accelerating and rotating 
charged black holes having also a NUT-charge has been presented. However, an accelerating 
NUT solution without rotation has not been identified yet within that class JZj . It is the goal 
of this section to try to construct such an accelerating solution. We shall consider next the 
uncharged C-metric, respectively the accelerating Zipoy-Voorhees solution ^ as the seeds 
in our solution-generating procedure. 

Expressed in the form given in the C-metric takes the simple form: 



2 , dy ^ dx 



where F{^) = 1 + 2mA^. We restrict our attention to case in which < 2mA < 1 and, in 
order to preserve the signature of the metric we restrict the values of the coordinates such 
that: 

-^<y<-i, -i<x<i. (21) 

In terms of these coordinates, spatial infinity corresponds to x = y = —1, the black hole 
horizon is located at y = —3^' "^^^^^ acceleration horizon corresponds to y = — 1. The part 
of the symmetry axis joining the black hole horizon with the acceleration horizon is x = 1, 
while the one joining the black hole horizon to infinity is x = — 1. 

In order to apply our solution generating technique we need to write the C-metric in 
Weyl form. Using the results from |16j we obtain: 

ds^ = -e-'f'dt' + e"^ [e'^idp^ + dz^) + p'd^''] , 

^-^ ^ {y'-i)ny) ^ ^2, ^ {y" - ^ (22) 

A^{x-yY ' f{x,y)G{x,yy 

where 

fix,y) = (y^-l)Fix) + il-x')Fiy), 

G{x,y) = [l + mA{x + yyY -m^A^il-xyy, (23) 
while the canonical Weyl coordinates p and z are defined such that: 

2 ^ {y^-l)il-x^)F{x)F{y) ^ (1 - xy)[l + mAjx + y)] 

^ A4(x-y)4 ' ^ A^{x-yy 

2 , r 2 fix,y)G{x,y) [ dy^ dx^ 



"^P ^"^^ A^{x-y)^ \iy^-l)F{y)^ {l-x^)F{x)' ' ^^^^ 
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Consider now the scaling transformation {ip,fi) — > (7'i/', where 7 is a real parameter. 
Applying it to the C-metric we obtain a new vacuum solution of the form: 



Applying it to the C-metric we obtain a new 



ds^ = - 



X 



■^^Fix) 



[f{x,y)G{x,y)f-' W ' ^Mv) 



dx"^ 



■dip' 



(l-x^)F(x 



(25) 



It is clear that by taking 7 = 1 we recover the initial C-metric. On the other hand, let us 
consider the zero-acceleration limit of this metric. Performing the coordinate transforma- 
tions: 



X 



cos 9, y 



1 



(26) 



while taking the limit A ^ and rescaling the metric by a constant factor A^'^"^ it is 
readily seen that we recover the Zipoy-Voorhees solution ([3j). Therefore, one could naively 
interpret the metric as describing an accelerating version of the Zipoy-Voorhees solution. 
However, the fact that the above metric is not the 'proper' accelerating Zipoy-Voorhees 
solution can also be seen from the fact that the 7 = 2 of this family should reduce to the so- 
called accelerating Darmois solution. This is the coincident limit of the accelerating Bonnor 
dihole solution that was recently found by Teo in j^j. In fact, a different metric describing 
the accelerated Zipoy-Voorhees solution, however, written in a very symmetric form has been 
presented by Teo in the same work. The 'proper' accelerating Zipoy-Voorhees solution reads: 



ds^ 



-e-^dt^ + [e2^(dp2 

{y'-l)F{y) f F{y) \ 
A^{x-yy \F{x)J 



+ dz^ 



+ p'^dip'] , 

{y' - l)F{y) F{yY 



(27) 



^^F{x 



fix,y)Gix,y) G{x,yY' 



where the canonical Weyl coordinates are again defined in (jzlj). Indeed, we see that the 
Darmois solution {i.e. the a = 2 member of this family) is clearly different from the 7 = 2 
member of ()25|1 and therefore we cannot actually interpret (j25|l as being an accelerating 
version of the Zipoy-Voorhees family.^ Nonetheless, since the metric fl25|l is only an interme- 
diate result in our solution-generating technique, we shall not further discuss its properties 
at this point, but limit ourselves to notice that one can apply the scaling transformation on 
Teo's solution and further generate a new family of vacuum metrics indexed by two distinct 
real parameters: 



^ + [e^^'idp' + dz') + p'dip'] , 



(2^ 



{y'-l)F{y) f F{y) \ 
A^{x-yy \F{x)J 



a-l 



,2^ 



{y'-l)F{y) F(y)"^"^F( x)("-^)^'' 
f{x,y)G{x,y) G{x,y) 



*We thank Edward Teo for pointing this out to us. 
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The next step is to charge the solution (PHj) using a general SL{2, R) transformation. 
Using the formulae (|12p we obtain: 



—e 



de + e^H^{x, y) [e^^'idp^ + dz^) + p^d^^] , 



A 



ACS 



(1) 



A2(a;-y)2 / F{x) 
{y'^-l)F(y) \F{y) 



a-1 



-dt, H^{x,y) 



^ " ' A'2{x-y)'2 \F(x) 



a-1 



7\ 2 
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Let us consider a few limiting cases of the above metric. Taking 6 Q and C — oo 
while keeping the product C^d constant we recover the uncharged metric ()28j) . On the 
other hand, the 7 = 1 member of this family should correspond to the charged accelerating 
Zipoy-Voorhees solution (see also |IH])- In particular, for a = 1 this should reduce to a 
charged version of the C-metric. However, unlike the known form of the electrically charged 
C-metric, in general, the above solution has a curvature singularity located at the roots of 
H^{x, y) = 0. Therefore its interpretation as a new form of the charged C-metric is dubious. 

By dimensionally reducing this solution down to three dimensions and dualising the 
scalar field x to an electromagnetic field as described in section |21 we find that the magnetic 
potential is given by: 



A., 



7 (1 — x^) {aF{x) + (1 — a)F{y)) 2m^ax 



C 



A^{x - y f 



AC 



(29) 



while the metric remains unchanged in this process. Finally, taking 7 = 2 and using ()20j) we 
find: 



ds'' 



{y'-l)F{x) fF{y) 



A^{x-yf \F{x 
, (F(x)F(y))"("-^) 



G{x,yy 
where we defined: 



C^ 
H{x,y) 

dy' 



{dt + A^d^f + 



H{x,y) 
A^{x - yy 



\F{y)J 



iy' - l)Fiy) 



+ 



dx"^ 



'l-x'^)Fix] 



x')F{y)d^' 



(30) 



H{x,y) 



1-6 



{y^-l)F{x) F{y) 



A^{x-yy \F{x) 



(31) 



This is the main result of this section. In the limit 5 0, C — 00 (with C'^5 constant 
and rescaling the t coordinate by a constant factor) we recover the Euclidean form of the 
accelerating Zipoy-Voorhees solution (j28p . It is manifest that the C-metric corresponds to 
the a = 1 member of this family. Another interesting limit to consider is the zero-acceleration 
limit. In this case it turns out that performing the coordinate transformations: 



X 



cos 6*, 



y 



1 



5A^ 



(32) 
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we obtain a general family of vacuum Euclidean solutions, indexed by a real parameter a, 
which interpret as the Zipoy-Voorhees generalisation of the Eguchi-Hanson soliton: 



2a r- 



4^25(1-2^)" 



r(r — 2m) 



{dz 



r — m)2 — m2 cos2 9 
2ma 



{dr^ + r(r — 2m)d0'^) + r(r — 2m) sin^ 



C 



cos 6'(iv3) . 



(33) 



Indeed, if we take a = 1 and properly rescale the z coordinate to absorb some constant 
factor we obtain a spherically-symmetric metric: 



t{t — 2?7T/ ) 

ds^ = 6 „ \, —T:^{,dz + 4m cos Odipf + (r^ - (5(r - 2mf) 



dr"^ 



r{r — 2m) 



+ dn'^ 



We distinguish now two possibilities. If we set directly 5 = 1 we can cast the metric in the 
following form: 



ds' 



1 



i?4 ; 



[dz + cosOdifY 



dR^ 



R' 



4 + . 



dQ' 



after redefining i?^ = 4m (r — m) and cr = 4m. This is the well-known Eguchi-Hanson 
soliton TS". On the other hand, if 5 7^ 1 then by redefining the radial coordinate such that 
R'^ — = r"^ — 6{r — 2m)^ with (1 — 5)n'^ = 4m^5 and rescaling z we obtain: 



ds' 



i?2 



{dz + 2n cos OdipY + 



R^ 



n 



{r + ^) [r + uVs] 



dR^ + {R^ - n')dQ? 



which we recognise as the Euclidean version of the Taub-NUT metric. It is now possible to 
set again 5 = 1 and recover the extremal Taub-NUT solution. 

Another case of interest is the one that corresponds to negative values for 5. Setting 
6 — *• —5, from the general form of the metric dSHI), we obtain a metric with Lorentzian 
signature:^ 

l + 5(l-^)'"r/ rfr-2m) 



ds 



1 + 2^)2° 



(r — m)2 — m2 cos2 9 



{dt+ ——cosedipY. 



(dr^ + r{r — 2m)d9'^) + r(r — 2m) sin^ 9d^p^ 



(34) 



where now 6 takes positive values only. Consider now the a = 1 member of this family of 



solutions. After redefining the radial coordinate such that R + 



r 



6{r — 2mY , where 



^Notice that C is not an essential parameter and it can be absorbed by a constant rescaling of the z 
coordinate. 

^We also changed the notation z ^ t. 
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{1 + 5)N'^ = ASm?, we obtain: 

ds^ = "Zi ^ ^-{dt + 2Ncosdd^f + ^^^^ ^dR' + {R' + N')d^\ 

^ +^ [R±f^{RTN^/5) 

i.e. the Taub-NUT solution |2Uj with mass M = ±^^^^ and NUT charge A^. On the other 
hand, setting 6^0 and taking the hmit C ^ oo while keeping the product constant, 
the metric ()34|) is readily seen to reduce to the Zipoy-Voorhees metric. 

Therefore, we expect that ()30|) describes the accelerating version of the family ()33|). 
Computing some of the curvature invariants for this metric one finds that generically there 
is a curvature singularity located at the roots of H{x, y) = as long a.s x y. However if we 
consider negative values of 6 {i.e. replace 6 —5) in the above metric we obtain a vacuum 
solution with Lorentzian signature and furthermore, we find that H{x, y) > always (for 
X ^ y). 



4 Properties of the accelerating Taub-NUT solution 

1 member of this Lorentzian 



In what follows we will concentrate our attention on the a 
family. The metric becomes: 

{y^-l)F(y) 1 



ds' 



A^{x — yY H{x,y) 



dt + - 



H{x,y) 
'A^{x-yy 



[1 - x^)F{x)d<^^ + 



A'^{x — yY 
dy' 



+ 



+ 



2mx\ 



A 



■) 

dx^ 



dip 



(y2 - l)F{y) (1 - x^)F{x) 



(35) 



where we denote: 



1+5 



{y''-l)F{y) 
A'^{x—y)'^ 



H{x,y) 

Let us first notice that the C-metric, respectively the Taub-NUT metric are included as 
limiting cases in the above solution. Indeed, taking 5—^0 and C oo while keeping the 
product constant, after rescaling the time coordinate with a constant factor we obtain 
the uncharged C-metric solution. On the other hand, the zero-acceleration limit is taken by 
performing the coordinate redefinitions and scalings of the parameters: 

x = cos6', y = s A'^6, t-^^, (36) 

in the limit A 0. It is readily seen that in this limit we obtain the Taub-NUT metric. 

To understand the properties of this solution it turns out to be more convenient to 
consider the above metric in Weyl form: 

ds^ = -e-'f'idt + A^d^f + [e^>'{dp^ + dz^) + p^d^p^] , 
'iCH .„ {y^ - l)F{y) 



,2m 



1 + 5e-2'^ ' 



f{x,y)G{x,y)' 



(37) 
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where the canonical Weyl coordinates p and z are given by ()24p and the expressions of e~'^ 
and e^^ in terms of the canonical Weyl coordinates are given in Appendix A. The general 
analysis of the above metric can now be done in parallel with the one corresponding to the 
uncharged C-metric. In particular, we see that if we restrict the values of m and A such that 
< 2mA < 1, then, in order to preserve the correct signature of the metric, the coordinates 
(x, y) have to take the range: 

-l<x<l. -^^<V<-1. (38) 

Now, it is well known that, in Weyl cylindrical coordinates, black hole horizons correspond 
to rods on the symmetry axis jlTj. Our interpretation of the above solution as describing 
an accelerating NUT-charged black hole will rely on the identification of such rods on the 
symmetry axis. 

We will define the symmetry axis to correspond to p = 0, i.e. it is the z-axis. Using (j^^ 
one sees that it corresponds to the four intervals: x = —1, x = 1, y = and y = —1. As 

we shall prove bellow, y = corresponds to the event horizon of the black hole, y = —1 

is the acceleration horizon, the line a; = 1 is the part of the symmetry axis between the event 
horizon and the acceleration horizon, while x = — 1 is the part of the symmetry axis joining 
up the event horizon with asymptotic infinity. 

To this end, notice that the asymptotic region x = y = —1 corresponds to z = ±00, 
while the end-points of the range of the coordinates (x, y) are mapped into z{x, y) as follows: 

/ 1 \ m f 1 \ m 1 

= ' ['^^ -2^)=-A^ '' = '[^^ - W )=T = = 2A^- 

Note now that e~'^\p=o vanishes at all the above points Zi, i = 1..3, it is positive for z < zi 
and Z2 < z < Z3 whereas both e~'^\p=o and e^''|p=o are zero for zi < z < Z2 and z > z^. 
We may then follow a similar analysis with the one performed in ^21j| to conclude that the 
regions Zi < z < Z2 and z > z^ are the Killing horizons of our accelerating solution. One 
can also see this by noting that the location of the horizons is given by the equation g^^ = 0, 
which in our case corresponds to the equation (y^ — l)F{y) = 0. Furthermore, by computing 
the area of each of the above horizons one can check that zi < z < Z2 has finite area and it 
corresponds then to a black hole horizon, while z > Z3 has infinite area and it corresponds 
to an accelerating horizon. Indeed, using the C-metric coordinates the area of the black hole 
horizon is readily found to be: 

2' "irm^ 



while the area corresponding to y = —1 diverges. 

Having determined that the above solution describes an accelerating object, let us turn 
now to a consideration of the 'cause' of the acceleration. The analysis of the conical singu- 
larities proceeds exactly as in the case of the uncharged C-metric. In particular, if we denote 
the periodicity of ip as Aip, then along a portion of the axis where the metric function e~'^ 
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is positive, e > 0, the deficit of conical angle will be given by^: 

A = 27r - Av9e"^|p=o. 



(40) 



Recall that if A < one has an excess of conical angle and this corresponds to a strut, if 
A > one has a deficit of conical angle that corresponds to a string, whereas if A = there 
is no conical singularity on that part of the symmetry axis. Since in our case the function e^'^ 
is precisely the same with the one corresponding to the uncharged C-metric, we deduce that 
in general there is a conical singularity residing in this solution and that, for appropriate 
values of Ayj, it can be chosen to lie along Z2 < z < {i.e. x = 1) or z < zi {i.e. x = — 1). 
In particular we find: 



One can remove the conical singularity on the segment Z2 < z < {i.e. x = 1) if one 
chooses Aif = 2tt/{1 + 2mA) but then there will be a positive deficit angle for z < zi 
{i.e. X = —1) and this can be interpreted as an semi-infinite cosmic string pulling on the 
black hole. Alternatively, for Aip = 2n/{l — 2mA) one can eliminate the conical angle for 
z < zi {i.e. X = —1) but then there will be a excess of conical angle for Z2 < z < z^ {i.e. 
X = 1). This is interpreted as a strut pushing on the black hole. The strut continues past 
the acceleration horizon and connects with the mirror black hole on the other side of it.^ 

Following the discussion in [12] let us consider next the presence of torsion singularities. 
In general these appear when the conical singularities possess a non-zero angular velocity, 
signified by a non- vanishing u = gtip/gtt along the symmetry axis. As is apparent from metric 
written in Weyl-Papapetrou form, near the symmetry axis p — > for a non-zero value of uj 
the coordinate will become a timelike coordinate and this will lead to the apparition of 
CTCs sufficiently close to the axis. In general, these CTCs can be eliminated only when u 
takes the same constant value along the entire axis of symmetry as in that case it is possible 
to perform a global coordinate transformation t — > t — uj\p=o(p to give a metric without such 
pathologies. For our accelerating NUT solution we find that on the symmetry axis p = 0: 



and therefore at the first sight there are unavoidable torsion singularities associated with 
this metric. However, one can still perform a coordinate definition = t+ on the line 
X = 1 respectively ts = t — near x = —1. Since the coordinate ip is periodic, this will 
introduce a periodicity for the time coordinate. However, this is precisely what is expected 
in the case of a NUT-charged solution [22] ■ 

5 Conclusions 

In this paper we constructed new solutions of the vacuum Einstein field equations in four 
dimensions via a solution generating method utilizing the SL(2, R) symmetry of the dimen- 

^The measurement of the proper circumference and proper radius must be performed in a frame for which 
the proper time dr = dt + A^dip = 0. 

^The existence of the second black hole on the other side of the acceleration horizon is obscured by the 
use of the Weyl coordinates. 




(41) 



Uj\x=±l — ± 



2m 
AC' 



(42) 
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sionally reduced action in three dimensions. Our method was based on the simple observation 
that a static axisymmetric metric as written in Weyl-Papapetrou exhibits a simple 'scaling' 
symmetry that allows one to generate a family of new static vacuum axisymmetric solu- 
tions, indexed by a real parameter. In particular, using this scaling symmetry one can easily 
generate the Zipoy-Voorhees solution from the Schwarzschild solution. 

We also made use of a charging method for static vacuum metrics, which dates back 
to Weyl p. We demonstrated a simpler alternative derivation of this transformation by 
using a SL{2, R) symmetry of the reduced Lagrangian in three dimensions. However, unlike 
previous applications of this transformation, we showed that with our simplified mapping 
and by combining this charging method with the scaling property, one is able to generate new 
vacuum stationary axisymmetric metrics. The Lorentzian version of the generated solutions 
gives a Zipoy-Voorhees like generalisation of the accelerating Taub-NUT solutions, while 
the Euclidean version gives a non-trivial two-parameter generalisation of the Eguchi-Hanson 
solitons in four dimensions. We focused our attention on a particular member of this family 
and we showed that it describes the accelerated version of the Taub-NUT space. 

As avenues for further research, it would be interesting to study in more detail the 
connection of the singular charged C-metric that we obtained with the usual form of the C- 
metric. In particular, it would be interesting to find a proper dilatonic generalisation of the 
accelerated Zipoy-Voorhees metric, one that would reduce to the proper charged C-metric 
in the appropriate limit. 
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A: The Weyl form of the accelerating Taub-NUT metric 

Following Emparan and Reall [2] we introduce the notation: 



Q = z-z,, R, = ^p^ + Ch = + R,R, + QCr (A-1) 

It can then be shown [2j that: 

_ ^ {y^-l)F{x) ^ + Ci = ^^-'^^^P^y\ R,-(,= + 

p I . _ {x + l){y-l)F{y) _ + 

A-{x-yy ' ^^-^^ 

while: 

_ {x-l)iy-l){2mA-lf _ {x - l){y - l)F{x)F{y) _ 2F{x)F{y) 
2A\x-yY ' 2A\x-yf ' A\x - yf 
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Then the Weyl form of the uncharged C-metric corresponds to the following expressions: 

R2 — (2 ' 

2, ^ 1 YuY2siRi-Ci)iR3-C3) 

A{2mA-iyRiR2R3 ^13 R2 - C2 ' ^ ^ 

from which we can readily find e""^ in (jHTj) . Finally, expressing x in terms of p and z we find 
US]: 



, _ 1 f {Ri + Ci){R2 - C2) , 2m Fi + F2 ^ 



where: 



Fo = Am^ARi + m(l + 2mA)i?2 + m(l - 2mA)i?3, 
Fi = -4m/?i - 2m(l + 2mA)R2 + 2m(l - 2mA)i?3, 

F2 = ^(l-2m2A2). (A-5) 

This completes the derivation of the Weyl-Papapetrou form of the accelerated Taub-NUT 
solution. 
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